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Abstract 

The present work is devoted to the evolution of random solutions of 
the unforced Burgers and KPZ equations in d-dimensions in the limit of 
vanishing viscosity. We consider a cellular model and as initial condi- 
tion assign a value for the velocity potential chosen independently within 
each cell. We show that the asymptotic behavior of the turbulence at 
large times is determined by the tail of the initial potential probabil- 
ity distribution function. Three classes of initial distribution leading to 
self-similar evolution are identified: (a) distributions with a power-law 
tail, (b) compactly supported potential, (c) stretched exponential tails. 
In class (c) we find that the mean potential (mean height of the sur- 
face) increases logarithmically with time and the "turbulence energy" 
E(t) = (v 2 (x, t)) = ((Vi/)(x, t)) 2 ) (mean square gradient of the surface) 
decays as t _1 times a logarithmic correction. In classes (a) and (b) we 
find that the changes in the mean potential and energy have a power-law 
time dependence, namely E(t) tx t~ p where the index p lies in the interval 
2 > p > (2 — d)/2. In class (c) the roughness of the surface, measured by 
its mean-square gradient, may either decrease or increase with time. We 
discuss also the influence of finite viscosity and long range correlation on 
the late stage evolution of the Burgers turbulence. 

PACS number(s) : 47.27.Gs, 02.50.Ey 
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1 Introduction 



The multi-dimensional Burgers equation 



;)j v vVv)=i/Av, (1) 

is a generalization of the well-known Burgers equation 

dv dv d 2 v . . 

-dl +V d-x =V d^- (2) 

The nonlinear diffusion equation (^), which was originally introduced by 
J. M. Burgers in f939 as a model of hydrodynamical turbulence |^], was later 
shown to describe a variety of nonlinear wave phenomena arising in the theory 
of wave propagation, acoustics, plasma physics (see, e.g., ||, ||, ||]). 

The Burgers equation (|^) describes two principal effects inherent in any 
turbulence Q: the nonlinear redistribution of energy over the spectrum and 
the action of viscosity in small-scale regions. Although external forces are not 
present in (^), the one-dimensional Burgers equation does describe the decay of 
turbulence, i.e. the nonlinear transformation of the random initial field v$(x). 
The Burgers equation shares a number of properties with the Navier-Stokes 
equation, namely the same type of nonlinearity, of invariance groups, of the 
energy-dissipation relation, of the existence of a multidimensional version, etc. 
However, Burgers equation is integrable and therefore is not sensitively depen- 
dent on initial conditions. The differences between the Burgers and Navier- 
Stokes equations are as interesting as the similarities JjJ and this is also true for 
the multi-dimensional Burgers equation ([l]). 

The three-dimensional form of equation ([!]) has been used to model the 
formation of the large scale structure of the Universe when pressure is negligi- 
ble, that is, during the nonlinear stage of the gravitational instability arising 
from random initial perturbations ||, [| |l(], pT[ . Other problems leading to the 
multi-dimensional Burgers equation, or variants of it, include surface growth 
under sputter deposition and flame front motion [ fL2| . In such instances, the 
potential "0 corresponds to the shape of the front's surface, and the equation for 
the velocity potential ip is identical to the KPZ (Kardar, Parisi, Zhang) equa- 
tion [jl2| [l3|, [l4|. For the deposition problem v = — VV> is the gradient of the 
surface. The roughness of the surface, measured by its mean-square gradient 
E(t) — ((Vt/>( x 7 t)) 2 ) — ( v2 ( x ^))i ma Y either decrease or increase with time. 
Nevertheless we will use the expression "turbulence energy" for this value of 
E(t). With external random forces the Burgers and KPZ equations describe 
phenomena such as turbulence without pressure, disordered systems, directed 
polymers, etc @ |I§ |TJ, [It], [l§. 

Here we will consider the evolution of the velocity field v(r, t) and potential 
^(r, t) as given by the noiseless multi-dimensional Burgers and KPZ equations. 
In this case the evolution of the fields is fully determined by the statistical 
properties of the initial field vo(r), which is assumed to be potential vo(r) = 
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— V^o(r). Our main attention will be given to the case of vanishing viscosity 
(y — ► 0), when the dissipative effects are important only in the vicinity of shocks. 

It is known that the asymptotic behavior (t — > oo) of Burgers turbulence 
strongly depends on the behavior of the structure function of the initial poten- 
tial d^,(p) = ({ipo( r + p) — ^( r )) 2 } a t large distances B |1. If the structure 
function (p) increases as a power law in space then the initial potential field 
is Brownian, or fractional Brownian motion, and some scaling may be used 
j| ||, pl|, [lS^I . In this case the turbulence is self-similar. The evolution of the ex- 
ternal scale L(t) of the self-similar solution in time is determined by the spatial 
behavior of d^(p) at large distances. Very recently a complete solution of the 
one-dimensional Burgers equation, with initial Gaussian white noise distributed 
data in the inviscid limit, has been obtained pof . 

If the structure function of the potential is bounded at p — > oo, then scaling 
arguments no longer work. Kida has proposed using, for the one-dimensional 
Burgers turbulence, a discrete cellular model of the initial conditions with inde- 
pendent distribution of potential in different cells. He has shown that the energy 
decays like with a logarithmic correction when the probability distribution 
of the potential in each cell has a stretched exponential tail ~ exp(— H"). Later 
it was shown |^2|, |2^, |24| that for an initial continuous Gaussian field all of the 
statistical characteristics of one-dimensional Burgers turbulence become self- 
similar and the energy decays as i _1 /(lni)2. Several models of the evolution 
of Burgers turbulence with an initial perturbation of non Gaussian type have 
also been proposed (O] and references therein, |25|, [2(| Q). It has 

been shown that the law of energy decay strongly depends on the statistical 
properties of the initial field with homogeneous potential. 

It is known that in the limit of vanishing viscosity the solution of Burgers' 
equation is reduced to searching for the absolute maximum of some function 
of the initial potential pl| . Consequently, the statistical properties of Burg- 
ers' turbulence are determined by the behavior of the probability distribution 
of the initial potential. One of the important results of the classical theory of 
extrema is that there are only three universal classes of "extreme value dis- 
tributions" of the sequence of independent and identically distributed random 
variables |3^]. This result was used in 1 33 for dealing with the problem of the 
equilibrium of low-temperature physics of disordered systems and partly for the 
onc-dimcnsional decaying Burgers turbulence. One of our main tasks is to show 
that, in the discrete cellular model of the initial condition for the d-dimensional 
Burgers turbulence, there are also three classes of universal self-similar evolution 
of the velocity and potential fields. 

The paper is organized as follows. In Section 2 we formulate our problem 
and list some elementary results about the Burgers equation. We also consider 
the evolution of potential and velocity fields for a simple kind of perturbation in 
individual cells at the initial stage. In section 3 we derive the general expression 
for the energy and probability distribution functions of the potential and vector 
velocity fields. In Section 4 we show that three types of initial condition lead 
to asymptotic self-similar behavior of Burgers' turbulence. Section 5 presents 
concluding remarks. We also discuss the influence of finite viscosity and long 
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range correlation on the late stage evolution of Burgers' turbulence. 



2 Description of the basic model 

2.1 Basic equations and local self— similarity 

We will discuss the initial- value problem for the unforced Burgers equation (|l|) , 
and consider only the potential solution of this equation, namely 

v(x,i)=-W(x,i). (3) 
The velocity potential ^(x, t) satisfies the following nonlinear equation: 

dtii> = ^(Wf + wV. (4) 

This equation is the same as the KPZ equation [|l2[ [lj, which is usually 
written in terms of the variable h = The parameter A has the dimen- 

sions of length divided by time and is the local velocity of the surface growth. 
Henceforth h has the dimensionality of length and is the measure of shape of 
the surface. Using the Hopf-Cole transformation ip = In [7 |3l], one can 
reduce (|4|) to a linear diffusion equation. We are mainly interested here by so- 
lutions in the limit v — * 0. Use of Laplace's method then leads to the following 
"maximum representation" for the potential and velocity fields in the limit of 
vanishing viscosity [|, [llj : 

f) = max$(x,y,t) , (5) 
*(x,y,t)=^(y)-fc^, (6) 

v(x,t) = ^^. (7) 

Here ipa(y) is the initial potential and Vo(x) = — V'0o( x )- In (Q) y(x, t) is 
the Lagrangian coordinate where the function $(x, y,t) achieves its global or 
absolute maximum for a given coordinate x and time t. 

At large times the paraboloid peak in @ defines a much smoother func- 
tion than the initial potential ipo(y)- Consequently, the absolute maximum of 
$(x, y,t) coincides with one of the local maxima of ipo(y). The Lagrangian 
coordinate y(x, t) then becomes a discontinuous function of x, constant within 
a cell, but jumping at the boundaries g, [llj. In each cell fluid particles move 
away from the small region near the cell center y^. The velocity field v(x, t) 
has discontinuities (shocks) and the potential field ip(x.,t) has gradient discon- 
tinuities (cusps) at the cell boundaries; these shock surfaces or walls form a 
connected structure. Inside the cells the velocity and potential fields (surface 
shape) have a universal self-similar structure: 

v(x,<) = ^-^, (8) 
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^(x,t) = v»o(y fc )- v 2 ; ; • (Q) 

The longitudinal component of the velocity vector v(x, t) consists of a sequence 
of sawtooth pulses with random positions of the shocks and "zeros" , just as in 
one dimension. The transverse components consist of sequences of rectangu- 
lar pulses with random amplitudes and random positions of the shocks. Wall 
motion results in continuous change of cell shape with cells swallowing their 
neighbors and thereby inducing growth of the external scale L(t) of the Burgers 
turbulence. At large times the behavior of the turbulence will be determined by 
the statistical properties of the initial potential and, moreover, by the statistical 
properties of local maxima V'o(yfc)- 



2.2 The cellular model and initial stage of evolution 

We consider the cellular model of the initial conditions, in which we assume 
that space is divided into identical cells, each having a volume dV = Lq, where 
Lo is the length of the cell and d is the spatial dimension. We assign an initial 
value for the potential which is chosen independently within each cell. The same 
approach was used in plj in order to find the probability distributions of the 
amplitudes and velocities of shocks for the one-dimensional Burgers turbulence 
and in ^ || for multi-dimensional Burgers turbulence in the case where the 
initial potential had a a stretched exponential tail. In the "shot-noise model" 
| pH it is assumed that the initial potential is a sum of potentials of "unit non- 
homogeneities" with random amplitudes and scales and with a Poisson ensemble 
for the position. 

First we consider the evolution of the velocity and potential fields of the 
"unit non-homogeneities" inside the cells assuming that there is no interaction 
between cells. We assume for simplicity that inside each cell the initial potential 
and velocity are isotropic and 

V>o(x) =Ml-x7©, |x|<Z , (10) 
v (x) = Vqx/Iq , |x| <l ,Vo = 2^ /% , (11) 
where "0o( x ) — 0, vq(x) = when |x| > Iq and Vq is the amplitude of the velocity 



at the border |x| = Iq. From the solution (g|7|) we have for the velocity field 

V(x ' t)= ^o(l + X W |X| <Xs ' (12) 
where x s is the position of the shock surface 

x s = l (l + V t/l Q ) 1/2 = (Ig + 2t/V) 1/2 ■ (13) 

In the one-dimensional case this solution is a well known N-wave || . 
The energy of the velocity field is given by 

E{t) = J v 2 (x,t)d d x = J vV-'dr, (14) 

o 
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and from (12 T^) it follows that 

m-*v+wM+. ^= rW2)( ;; 2) . as) 

where T(z) is a gamma-function. 

From ([l5]) we see that for d = 1 the energy £?(£) decreases with time, for d = 2 
-E(i) is constant, and for d > 3 -B(i) increases with time. It is possible to show 
that E(t) may increase with time when d > 3 even if the viscosity coefficient ^ 
is finite. The increase of the energy in the multi-dimensional Burgers equation 
(in contrast with d = 1) is the result of this equation not having a conservation 
form. 

Moreover it easy to show that in the multidimensional Burgers equation 
the energy does not conserve in the limit of vanishing viscosity even at the 
initial stage of evolution, when the velocity field does not have shocks. Let 
us consider the evolution of the continuos velocity field having only a radial 
component u r (x, f) = v(r, i). In the limit of vanishing viscosity we have the 
implicit solution of Burgers equation v = vq(t — vt). Substitute this solution 
into ( [hi] ) and replace the Eulerian coordinate r by Lagrangian coordinate y: 
r = y + tva (y) , we have the following expression for the energy 

oo 

Eit) = ^1 J vtm + tvoiy))^ (l + *^) dy. (16) 
o 

From (|l6|) we see that the energy E(t) of continuos field conserves only in one 
dimensional case (d = 1). In particular for d — 2 we have from ( [l6| ) 

oo oo 

E(t) =2irj v 2 (y)dy + t j v 3 (y)dy . (17) 
o o 

Thus for the positive perturbation of the initial velocity the energy (mean- 
square gradient of the surface) increases with time, and for the negative pertur- 
bation of the velocity the energy decreases with time. This effect takes place 
due to the relatively more significant contribution into mean-square gradient of 
the surface the regions with large distances r (|l6|). 

The asymptotic behavior of the potential and velocity fields is determined 
only by the positions and values of the maximum of the initial potential in the 
cells (|||^). For relatively large time (Vot ^> Iq) the velocity inside the cells has 
the universal behavior 

v(x, t) = x/t , |x| <x s = (2^ t) 1/2 • (18) 
The position of the shock surface x s and the energy of the field 

±tl d-2 

E(t)~%* t— (19) 
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arc determined only by the value of the initial potential maximum ipQ. It is 
easy to show that this asymptotic behavior takes place for the arbitrary initial 
localized perturbations with the potential maximum ipQ. 

Let us assume (Lq 3> Iq) and that the maximum value in the "i-th" cell is 
ipo.i- The interaction between cells begins when the shock surfaces of neighbor- 
ing cells come in contact, i.e. when \x St i\ + \x s ,j\ = Lq. At rather large time the 
border between the cells becomes a hyper-plane, moving with constant velocity 
IV-'cm — V'Oj |/|yi — Yj I away from the cell with a larger value of the unital poten- 
tial ||. When the potential maximum in all of the cells is the same, i.e. when 
ipo,i = V'Oi then at large times the hyper-surfaces between neighboring cells are 
immobile and the velocity v(x, t) and potential ^(x, i) fields have stationary 
periodic structure. The gradient of the velocity field decreases as t -1 , which 
means that the energy density E(t) = (v 2 ) decreases according to 

E(t) ~ L 2 /t 2 . (20) 

For a random distribution of the potential maximum we have the per- 
manent absorption of cells with low values of the potential so that the external 
scale of the turbulence L(t) increases with time. The mean energy evolution 
strongly depends on the statistical properties of the maximum potential distri- 
bution. The energy decay is slower than for the periodic structure ((2^). The 
interaction between cells changes also the law (|l9|). Moreover, for d > 3, due 
to the interaction between cells the increase of energy E(t) ( pj| ) may change to 
a to a decrease. It will be shown later that asymptotically E(t) oc t~ p with p 
lyings in the interval 2 > p > (2 — d)/2. 



3 Energy evolution and probability distributions 
of the turbulence and the interface 

We assume that the initial random potential ipo (y) is statistically homogeneous. 
Then, from equation (||) it follows that the turbulence energy is determined by 
the time derivative of the mean potential 

E(t) = <v 2 (x,*)> = <(W,(x,t)) 2 > = 2^|^ . (21) 

Let (H, t) denote the probability density of ip at time t and let (H, t) 
denote the cumulative probability to have ip < H, given by 

H 

Qi,(H,t) = J P^(A,t)dA. (22) 

— OO 

The mean value of the potential (ijj) (i.e. mean value of the height of the surface) 
at time t is then expressible as 

W = J HP^{H,t)dH. (23) 
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Thus we can find the mean energy of the turbulence if we know the probability 
distribution function P^(H,t), where P i f,(H,t)dH is the probability that the 
absolute maximum of ^(x, y, t) (^|) lies in the interval H,H + dH for all y. 

The function P^^ y (H, y; x, t) is the joint probability distribution function 
(p.d.f.) of ip and y and P^,dHd d y is the the probability that the absolute 
maximum of $(x, y, t) lies in the interval H,H + dH in the volume element 
d d y. By integrating the P^ iy (iJ, y; x, t) over H we obtain the p.d.f. of the 
Lagrangian coordinate y: 

+00 

P y (y;x,i)= J P^ y (H, r ,yL,t)dH , (24) 

—00 

which permits one to find the probability distribution function of the velocity 
field tff). 

Since the initial potential is homogeneous we will consider the point x = 
and no longer include x in the parameter of the probability distribution function. 
In the "cellular" model we assume that the space is divided up into elementary 
cells of volume dV — Lq. We consider the late stage of evolution when the 
external scale L(t) is much greater then Lq. On this scale the initial potential 
field is homogeneous. (Formally we can introduce a random uniform distribution 
of cell positions.) 

The asymptotic behavior of the probability distribution function of abso- 
lute maximum of a large number of random quantities is determined by the 
asymptotic properties of their cumulative distribution function at large values. 
We assume that in each cell the cumulative probability function of the initial 
potential ^o(y) is the same and can be represented in the following form: 

Prob{ip Q <H) = F(H) = 1 - f(H) . (25) 

The cumulative probability distribution function for $) (||) in the i-th cell is 

Fi(H) = 1 - fi(H) , f i (H) = f(H+?±). (26) 

Here yj is the coordinate of the i-th cell. Due to the independence of the initial 
potential in different cells we can determine the cumulative distribution function 

ip as 

Pro6(V < H) = Q(H, t) =\\F 1 {H) 

= -fi(H)). (27) 

From ( p6|) we see that fi(H) decreases with the distance |y|. 

We will consider the limit when t becomes large, so that the number of 
factors in ( p7| ) becomes significant. This means that the absolute maximum 
of $ is reached for rather large H, so in each cell we have fi(H) << 1, thus 



1 — fi(H) ~ exp [-fi(H)], and we can rewrite (27) in the form 

Q(H,t) =e- N ^, (28) 
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where 

N(H, t)=J2 fi{H) =T,f( H + Yt ] - (29) 

i i 

At large t when the difference between /j in neighboring cells is small, we can 
replace the summation in (^9[) by an a volume integral: 



N{H,t) = jjJ f(H+^)d d y. 



(30) 



For large H the events for which tpi is greater than H are rare and they have a 
Poisson distribution. Thus N in (|2^) is a mean number of events when ip > H. 
From (E8Lp0ft we have two equivalent expressions for the mean potential: 



+ oo 



m= J H-^e- N ^dH= J H(e,t)e- e dQ, (31) 

-oo 

where H(Q, t) is a solution of the equation 

N(H,t) = Q. (32) 

Formally we can integrate in ( |3l"| ) over the infinite interval of H since the prob- 
ability of H being negative is negligibly small when t is large. 

In order to find the probability distribution function of the velocity we first 
of all need to find the joint probability distribution of value and the coordinate 
of the absolute maximum of $. The probability that the absolute maximum is 
in the i-th cell, with a value in the interval H,H + dH, is 



Prob(y eAj,He[H,H + dH]) 
d 

~ ~ dH 



■^fAH)Y[(l-h{H)). (33) 



Here the product over all cells with i ^ j is the probability that outside the j-th 
cell all local maxima are less than H. 

With the same assumption as before we can find the joint probability distri- 
bution coordinate and value of the absolute maximum. At large t we have 

P y (y,t) = P(y&A j )/L d . (34) 

and after integration over H we have the following expression for the probability 
distribution function of the absolute maximum coordinate 
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Integration by parts leads to 

+00 

— 00 

= zgM* (e)+ ») e_ed8 - (36) 



Here H(Q) is a solution of equation (|32|). Using ( |30| ) it is easy to see that the 
norm of P y (y,t) is unity. 

In the theory of the absolute maximum of random sequences p2[ the large 
parameter of the theory is M - the number of "points" in a sequence. In our case 
we consider formally an infinite number of points, but owing to the parabolic 
form used in (|2^) the effective number of cells is finite and increases with time. 
We stress again that expressions ( |30| , |36| ) are valid only when L(t) » Lq. 



4 Three types of self-similarity for Burgers Tur- 
bulence 

In a previous section we derived expressions for the probability distribution 
function of the potential ^>(x, t) (^,^) and for the mean potential (|3l| ) which 
determines the energy evolution of the turbulence (pl|). We also derived the 
form of the probability distribution function for the absolute maximum coordi- 
nate ([35j) which coincides with the form of the velocity probability distribution 
function ((?]). It is known that there are three types, known as the Frechet, 
Weibull and Cumbel classes, of asymptotic behavior for the probability distri- 
bution function of the large sequence of random quantities [^|[ |33|. In this 
section we will show that these three types of initial condition lead to three 
types of asymptotic self-similar behavior of the Burgers turbulence. 



4.1 Distributions with a power- law tail: Frechet class 



Let us now assume that the cumulative probability distribution function of the 
initial potential (25) has a power law tail 



/ (//)= 7 



(37) 



In the integral ( pOj ) for N(H,t) we can integrate over the radius r = |y| and 
obtain 



r(d/2) 



r.d-1 



dr 



2Ht 



(38) 
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This integral converges when 7 > i , and gives the asymptotic distribution of 
the potential ip when the probability distribution function of the initial potential 
decays sufficiently rapidly. Replacing the variable r by r — (2tH) 1 ^ 2 x allows us 
to write N(H, t) as a power-law function of H 

H 1 t d / 2 



where Cd depends on d as well as 7. Thus we have the self-similar cumulative 
function ( p8| ) of the potential (height of the surface) is given by 



Q(H) = exp 



H \t 



(40) 



Here we introduce t n i as a nonlinear time of evolution 

*<* = inm ■ ( 41 ) 

It follows from (^) that the probability distribution function of potential is 
self-similar with the amplitude scale H s (t): 

/ t \d/(2 1 -d) 

H s {t)=HA — \ . (42) 

\ t n l J 

The time evolution of the velocity variance is determined by the time 
derivative of the mean potential (|3l|). It is easy to see from ( pitf ) that the 

expression for the mean potential converges if 7 > — + 1, and in this case the 

mean value {iji) ~ H s (t). For 7 < d the mean potential (ip) increases faster in 
time then the local height of the surface which grows proportionally tX. The 
superlinear increase of mean potential is possible while the KPZ equation (|J) 
describes the evolution of shape ip(x., t)) ~ h(x, t) = H (x, t) — tX of the surface 
and the normal growing of the surface is excluded a for this equation (l^, [l3[ . 
From ( ^l|) the velocity variance (turbulence energy) is given by 

We stress here that for this expression to be valid 7 must satisfy 7 > 7 cr = — + 1. 
From ( ^3|) we see that when 7 lies in the interval 

^ + !<7<d, (44) 
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the velocity variance increases with time. Thus the velocity variance increases 
only if the spatial dimension d satisfies d > 3. For d = 3 we have the same 
results as the "shot-noise model" of the initial perturbation p5| . On the left 

hand side of the inequality (p4|), if 7 = — + s, (e < 1) we have for the index p 
in @ 

Pe = 2 • (45) 

This means that (a) when d = 1 the energy always decays with p ~ 1/2, (b) 
when d = 2 the variance decays very slowly with p = e, (c) when d = 3 the 
variance increases with p ~ —1/2, and (d) when d ^> 1 then p ~ (2 — <i)/2. 

This behavior of critical index has a simple explanation if we consider the 
velocity field evolution in an elementary cell at the stage before the interaction 
with other cells takes place (Section |2.2| ). It follows from dl9) that the energy 

E(t) of an elementary volume varies as E(t) ~ ipQ 2 t 2 . When the power index 
7 of potential probability distribution function is close to the critical value j c , 
then probability distribution function of the initial potential has a relatively slow 
decay and the general evolution of the turbulence variance is determined by the 
cell with very high amplitude. It is easy to see that at the critical point the 
time dependence of the mean variance ( f45| ) is the same as the time dependence 
of the energy of an initial cell (fl9|). 

As 7 increases the interaction between the elementary cells begins to play 
a more and more important role. When 7 = d and {d > 2) the increase of 
variance E(t) will switch to a decrease with time. When 7 > d the variance 
will decay as i -1 which is cl ose t o the law for the stretched exponential type of 



initial potential (see Section 4.2) 



It is easy to show from (|36| ) that the probability distribution function of the 
coordinate of the absolute maximum of the potential is also self-similar, i.e. 

w ® - ik Ws im> • (46) 

where W s (z) is the dimensionless form of probability distribution function 

00 

W s {z) = A 







1 2 



(47) 



A is normalization coefficient depending on d and 7. In (46) L(t) is the external 
scale of turbulence and 

t 2-,-d \ tnl / 

The distribution ( p7| ) decays ~ z 2 at small distances and by the power law 
at large distances: W(z) ~ A/z 2j . From this we see again that the mean energy 
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is finite when 7 > r y cr = — + 1 , which is also the condition that the mean energy 

in separate cells ( ^9|) E(t) ~ J H~z~ f'(H)dH be finite. Note that the mean 
initial energy © B(0) = (E }/L$ ~ J H 2 f'(H)dH is finite if 7 > 2 for all 
value of d. 

4.2 Compactly supported potential: Weibull class 

Let us now assume that the initial potential is bounded by some value H m 
(■00 (y) < H m in all cells) and that the cumulative probability distribution func- 
tion has a power law behavior in the neighborhood of the maximum, i.e. 

f(H)=(^^-Y, H<H M (49) 



and is zero for H > Hm. By the same procedure as above (see Section 4.1 ) 
it follows that the cumulative probability distribution function of potential has 
the form (128) with 



LqH™ 

o±tl A/2 ' 
2 2 Ti a l £ 





and is self-similar for arbitrary dimension d and power index a. The mean 
potential tends to the maximum Hm according to the power law 

(iP) = H M -cH*(t nl /t)*fc. (51) 

Here c is a positive numerical coefficient, and t n i is the nonlinear time t n i ([ill), 
where the constant Cd is determined by the expression (|50|). 

From (|5l]) it follows that the energy always decreases according to the power 

It follows from (|36|) , (fl9"|) , (|50|) that the p.d.f. of the absolute maximum coor- 
dinate y is self-similar with the spatial scale L(t) 

L(t) ~L {t/t nl )^FZ. (53) 

The general result for the length scale evolution for systems in the Weibull class 
has actually been noted in |29[| (although not explicitly derived). It is easy to 
see that the energy (^ , (|52j)~and the external scale (f4g|) , (|53"1) of the turbulence 
are described by the same expressions for a distribution having a power-law tail 
and for a compactly supported distribution if we set a = —7 in ( f49| ) (7 < 0). 
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The form of the distribution of the absolute maximum coordinate is now 
defined by the integral 



W s {z) = A J (e^ - z 2 y e- e d8 , 

e,(z) 

e4z) = z 2a+d . (54) 
It can be seen that this distribution decays exponentially at large z 

W s (z) ~ e~ z2a+d z a ( 2 - d - 2a K (55) 

It is interesting to compare the evolution of Burgers turbulence for a com- 
pactly supported initial potential with the evolution of turbulence for a Gaus- 
sian initial potential having scaling properties with dimension h, such that 
> h > — 1 For values of h that satisfy this condition the initial potential 
correlation function has a singularity at the origin. In a special class of ran- 
dom solution of the one-dimensional Burgers equation was constructed. For this 
solution the self-similar behavior is true at all times and not just asymptotically. 
Such scaling behavior corresponds to the characteristic length L(t) ~ i 1 /( 2 ~' 1 ). 
Comparing this result with ( |53| ) we see that for a = —1/h and a > 1 the evo- 
lution of turbulence with a compactly supported initial potential tends asymp- 



totically to the self-similar solution described in 30 . 

From ( |5^ ) we see that when probability distribution function of initial po- 
tential decays very rapidly (a — > oo) the mean energy decays as i _1 just as it 
does for the stretched exponential type of initial potential. 

The special case where a = 1 is equivalent to the uniform behavior of the 
probability distribution function for the initial potential over H < Hm, that is 

to say: W(H) = - 9 QJP = H' 1 . In this case the energy E(t) ~ i2(i+rf)/(2+d) 
oH 

which is equivalent to the asymptotic law for the uniform distribution of the 
initial potential |2£| . 

The limit a —* when the energy decays as t~ 2 in all dimension d is also 
interesting. In this case the probability distribution function of the coordinate 
(and of the velocity) has the universal form W s (z) ~ exp (— z d ) for all z. It easy 
to see that for a — > the integral scale of the turbulence does not increase with 
time. The case of a — » has a simple explanation if we rewrite the cumulative 
probability distribution function for the initial potential corresponding in the 
form 

f(H) = PE(H M - H) , 
W(H) = -?^p-=P6(H-H M ). (56) 

Here E{z) is the unit function, S(z) is the delta-function, and P is the proba- 
bility that in some cells the potential amplitude is exactly equal to Hm- This 
probability may by rewritten as P = (Lo/Lp) d where Lp is the typical distance 



14 



between cells when the initial potential is equal to -ff* and Lq is the size of 
the cells. It is easy to see that in this case the universal behavior of energy 
decay E(t) ~ t~ 2 is due to the fact that after some intermediate stage we have 
"frozen" turbulence. This means that the spatial structure of the velocity will 
not change in time and only the amplitude of the velocity will decay like AL/t, 
where AL is the random distance between two cells with equal initial potential, 
i.e. (Hi = Hj = Hm)- In the case d = 1 the velocity will be a sequence of the 
lines v = (x — Xi)/t with the immobile shocks at (2^+1 — Xi)/2 1 where Xi is the 
position of the point where tpi — Hm- This case was considered in [27j but 
with some other assumptions and with other tools. In arbitrary dimension d 
the energy decays as L 2 P /t 2 , and the probability distribution function of the La- 
grangian coordinate and velocity will have the stationary form W(z) = Ae~ zd 
with the scale L ~ Lp ~ LoP^ 1 ^. For d = 2 this distribution has a Gaussian 
form. In conclusion we note that this asymptotic "frozen" behavior for bounded 
initial perturbation of potential takes place for arbitrary a when d — > 00. 



4.3 Stretched exponential tail 

of the initial potential: Gumbel class 

Let us now assume that the cumulative probability distribution function of the 
initial potential ( |25| ) has a stretched exponential tail 

f(H) - (H/H p ) a exp l-(H/H*f] . (57) 

For initial conditions of this type of the distribution potential is localized in a 
narrow region AH near the mean value Hq ^ iJ*. In the integral (|3^) over y 
we can take into account only the quadratic term in the exponent and obtain 
the following expression 



2-KtH? 



-| d/2 



e 



PH^L 2 

Let us introduce the dimensionless potential z so that 



(58) 



ff = HJ T + oJ' (59) 

where ho is a solution of the transcendental equation N(H*ho) = 1: 

2 n + d(l-/3) / t \ d/2 h & 



Here t n i is the nonlinear time 



In— ) (60) 



tni -2^m{H: ■ (61) 



15 



The variable z in ( p9[ ) has a universal double-exponential distribution at t 3> t n ; 

F(*) = e~ e ~ z . (62) 

When t 3> the distribution of the potential tp is concentrated in the narrow 
region AH /H ~ h~ p ~ [dlnt/t,^ 1 near the mean value (H) = H^Hq ~ 
(In t) 1 /' 3 . From (|2l]) we have the following expression for the energy of turbulence 

From ( |63|) we see that the energy decays according to the universal law E(t) ~ 
t -1 with some logarithmic correction: the decay is faster if p > 1 (tail of 
probability distribution function f(H) decays faster then the exponential law), 
and slower if ft < 1. It must be stressed also that the law of decay does not 
depend on the dimension of the space and so is the same as in one-dimension 

IB- 

The probability distribution function of absolute maximum coordinate y and 
the velocity v (^) are Gaussian and the variance of each component is given by 



m 3 — 1 P- 1 ' 

K Qat/tra)-*- 

(v 2 )=L 2 /t 2 . (64) 

The two-point probability distributions of the velocity and correlation func- 
tions were found in || |8| for the special case p = 2. The shape of the two-point 
probability distribution function of the longitudinal component is the same in 
the space of different dimensions, and coincides with the analogous p.d.f. of 
the one-dimensional Burgers turbulence. The transverse velocity components, 
unlike the longitudinal ones, are statistically independent in different cells and 
have a Gaussian probability distribution inside them. These results may be 
extrapolated to the arbitrary stretched exponential tail of the initial potential. 



5 Concluding remarks 

The present work has considered the evolution of random solutions of the un- 
forced Burgers and equations in d-dimensions in the limit of vanishing viscosity. 
The main statistical assumption is the independence of the initial velocity po- 
tential in different cells. 

We show that the asymptotic behavior of the turbulence at large times is 
determined by the tail of cumulative initial potential probability distribution 
function F(H) = 1 — f(H). We show that three classes of initial distribution 
lead to the self-similar evolution of the turbulence at large times. In the theory 
of extremes these limiting distributions are known as (a) Frechet class when 
f(H) oc H-~i; (b) Wcibull class when f(H) cx (H max - #)~ 7 ,7 < 0; and (c) 
Gumbel class when f(H) oc expl—H 13 ). One can find in |32| more general 
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conditions which are necessary and sufficient for the probability distribution 
function f(H) to belong to each of three types. 

We show that the mean potential (mean height of the surface) increases 
with time in cases (a),(b) according to the power law (ip) oc At r , r — d/ (2-f — d) 
(A < for 7 < 0) while in case (c) it increases logarithmically according to 
(ip) oc (lnt) 1 /' 3 . For the Gumbel class the distribution of potential is localized in 
a narrow region near the mean value. For the Weibull class a p.d.f. of potential 
exists for 7 > d/2 and the mean value of ip is finite for 7 > d/2 + 1. The 
mean square gradient of the surface (turbulence energy) E(t) — ((Vip(x, i)) 2 ) = 
(v 2 (x, t)) has the power-law dependence E(t) oc t~ p ,p = 2(7 — d)/(2j — d) 
in cases (a) and (b) and decays according to i _1 (lni)( 1 ~' 5 )/' 3 in case (c). For 
relatively slow decay of the initial probability distribution of potential 1 + d/2 < 
7 < d and dimension of the space d > 2 the energy E(t) increases with time. 

We show that the p.d.f. of the velocity is self-similar with the scale L(i)/t, 
where L(t) is the external scale of the turbulence. In cases (a) and (b) the 
external scale of turbulence increases according to the pure power law L(t) oc 
t m ,m = 7/(27 — d) while in case (c) it increases like L(t) oc f 1 / 2 (m£)^ -1 )/' 3 . 
One can see that for fast decaying initial distributions of potential (I7I — > 00) in 
classes (a) and (b)) the law of external scale increase and the law of energy decay 
tends to the corresponding laws of the Gumbel class. In the special case 7^—0 
(class b) we have the "frozen" turbulence which means that the structures of 
the potential and velocity fields conserve L(t) = const and the amplitude of the 
velocity increases as t" 1 . We note that for the compactly supported potential 
this "frozen" behavior takes place for arbitrary 7 when d — > 00. 

We now discuss what influence finite viscosity has on the asymptotic behavior 
of the Burgers turbulence at large times. For large initial Reynolds number 
(Reo ^> 1) we still have the cellular structure of the turbulence at relatively 
large times. For finite Reynolds number it will be characterized by two scales, 
namely the external scale L(t) and inner scale 6 the latter being a typical width 
of the shock surface S . Owing to viscosity the inner scale increases as S ~ vt/L. 
However due to the increase of L(t) the relative width of the shock 8(t)/L(t) ~ 
vt/L 2 {t) ~ i?e _1 (i) may either decrease or increase with time. Here we write 
the Reynolds number as Re(t) = V(t)L(t)/v = L 2 (t)/vt, since the local slope 
of the velocity is 1/t and the maximum velocity is of order L(t)/t. 

It is easy to see that for classes (a) and (b) we have a power law for the 
Reynolds number, namely Re(t) oc t r ,r = d/(2 / y — d). Thus for the Frechet 
class (7 > d/2) the Reynolds number increases with time and hence even when 
the viscosity is finite we have the strong nonlinear stage of evolution at large 
times. For the Weibull class (7 < 0) the Reynolds number decreases with 
time and at large times the evolution of the turbulence will be determined 
only by the linear diffusion. This is shown in p8[ for the case of a uniform 



distribution of the initial potential. In 28 it was stressed that although one 
finds a universal power law growth for the energy decay (for a system with a 
bounded and flat initial distribution), the velocity- velocity correlation function 
has a more complicated form, and there is no simple dynamical scaling in the 
system in the 'non-linear regime'. For the Gumbel class (f(H) oc exp(— H")) 



17 



the Reynolds number Re(t) oc (lni)^ 1- ' 3 '/' 5 . Thus we may expect that for (3 < 1 
we have conservation of the cellular structure at large times. For (3 > 1 the 
nonlinear evolution is only an intermediate asymptotic which changes at large 
times to linear decay. This effect is considered in || in the case of a Gaussian 
initial perturbation. 

Let us now move to the case when the initial potential is long-range corre- 
lated or non-homogeneous. We consider first one-dimensional turbulence, as- 
suming that the initial velocity is homogeneous with a spectrum E v (k) ~ a 2 |fc|™ 
at small wavenumbers k and falling off quickly at large wavenumbers. For the 
spectrum of potential we have E^{k) ~ a 2 \k\ n ~ 2 . For a Gaussian velocity it is 
shown in J|5| that there are three regions of n with different behavior of turbu- 
lence. When — 1 < n < 1 the long-time evolution of the spectrum is self-similar 
and the external scale of turbulence increases as L(t) oc (at) m ,m — 2/(3 + n) 
and is determined by the "amplitude" of the large scale component a. When 
1 < n < 2 the spectrum has three scaling regions: first, a \k\ n region at very 
small k with a time-independent constant, associated with long-correlated re- 
gions in physical space, second, a k 2 region at intermediate wavenumbers which 
is related to the self-similarly evolving "inner region" in physical space and, 
finally, the usual k~ 2 region, associated to the shocks. The growth of the ex- 
ternal scale is now determined by two integrals of the initial spectrum and 
L(t) oc t 1 / 2 (lnt)~ 1 / 4 . Switching wavenumber from the \k\ n to the k 2 region 
tends to zero faster than the energy wavenumber ~ 1/L(t) and asymptoti- 
cally we have the self-similar evolution of the spectrum. For n > 2, long-time 
evolution is also self-similar and L{t) cx ^^(lnt) -1 ' 4 . Thus for a Gaussian 
perturbation we have one critical index n for the behavior of external scale and 
energy. When n — 1 the index of the power law dependence of L(t) continuously 
transforms from m = 2/(3 + n), (n < 1) to the m =1/2 index of the leading 
term in the region n > 1. 

Let us assume that for the multi-dimensional Burgers turbulence the initial 
potential is isotropic and has a power-law dependence at small wavenumbers, 
i.e. Ej,(k) ~ a 2 \k\ n ~ 2 . The variance of the potential is determined by (ijj 2 ) ~ 

oo 

J E x f ) {k)k d ~ 1 dk and is finite when n > 2 — d. For 2 — d < n < — d its structure 
o 

function is d^(x.) = ((^ ( x ) — V'o(O)) 2 ) ~ |x| 2-,i ~ d . Using the rescaling of the 
structure function in (||) we see that the external scale increases as L(t) ~ 
(at) m , m = 2/(2 + n + d) HfOJ- This law does not depend on the p.d.f. of the 
initial potential. Introducing the Reynolds number directly through the Hopf- 
Cole solution g Re(t) = L 2 {t)/v = V(t)L{t)/vt we see that in this region of n 
the Reynolds number increases with time according to oc £( 2 -"- d )/( 2 +™+ d ). it 
means that at late stage the turbulence has a strong nonlinear cellular structure. 

Considering the case of independent amplitudes in different cells, we see 
that the correlation function of the initial potentials equal to zero if the distance 
between the points is greater then the size of the cell. Thus, the discrete model of 
the initial conditions considered in the present paper is equal to the uncorrelated 
potential with n > 2. Here we show that the laws of evolution of the external 



18 



scale L(t) @,(||),(|64|) and the energy E{t) @,(§2|),(|6|) are very sensitive to 
the tail of the p.d.f. of the initial potential. For the stretched exponential tail 
(i.e. the Gumbel class) the leading term for the evolution of the external scale is 
t 1 / 2 . Thus we can expect that in multi-dimensional turbulence the long-range 
correlation of the potential for 2 — d < n < 2 does not play an important role 
for the evolution of the external scale and energy and we have only one critical 
index n. This is shown in |^4|, |35| for the one dimensional case with a Gaussian 
potential. When n = 2 — d the universal index of the power law dependence 
of L(t) continuously transforms from m = 2/(2 + n + d), (n < 2 — d) to the 
index to = 1/2 for the leading term in the region where n > 2 — d. In this case 
we have also the same critical index n = 2 — d for the energy. In the interval 
2 — d < n < 2 we have conservation of the velocity spectrum E v (k) oc \k\ n at 
very small wavenumbers. But this small region is not significant for the energy 
of the turbulence and asymptotically the spectrum of the velocity tends to the 
self-similar evolution. 

We have a much more nontrivial situation when the potential distribution 
has a with power-law tail (Frechet class) or when it is a compactly supported 
potential (Weibull class). For n < 2 — d we find, using the scaling properties of 
solution (||), that the external scale increases as L(t) ~ (at) m , to = 2/(2 + n + d). 
Thus the evolution of L(t) and the energy E(t) ~ L 2 (t)/t 2 do not depend on 
the p.d.f. of the initial potential. For the uncorrelated potential n > 2 the 
power indexes laws of the scale evolution m = 7/(27 — d) and energy evolution 
p = 2(7 — d)/(2-f — d) depend on 7. The case of a homogeneous continuous 
potential with a power index 2 — d<n<2is equal to the existence of long- 
range correlation of the the potential amplitudes in cells. We can expect that 
in the region 2 — d < n < 2 the long-range correlation of potential influences 
the energy decay. 
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